UNCLASSIFIED 


HepAtuiuced. 

^  Ute 


ARMED  SERVICES  TECHNICAL  INFWMAnON  AGENCY 
ARLINGTON  HALL  STATId 
ARLINGKW  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  govemment  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U»  S» 
Govemment  thereby  inc\irs  no  responsibility,  nor  any 
obligation  ;rtiat soever;  and  the  fact  that  the  Govem¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacttire,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


N 

a 


NtVWEPS  lEPIlI  7i4S 
lOTS  TP  2S63 
tOPY  75 

CALCULATION  OF  CRITICAL  TEMPERATURE 
AND  TIME-TO-EXPLOSION  FOR 
PROPELLANTS  AND  EXPLOSIVES 


By 

P.  A.  Longwell 
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ABSTRACT.  The  solutions  available  in  the  literature  for  the  critical 
temperature  of  various  geometric  shapes  and  boundary  conditions  are 
gathered  together,  and  a  nomograph  is  presented,  which  greatly  facil¬ 
itates  numerical  calculation  of  the  critical  temperature  for  specific 
cases. 

Nomographs  are  also  presented  for  the  numerical  calculation  of 
the  time  required  to  cause  explosion  in  the  event  the  temperature  at 
the  surface  of  an  object  is  held  at  a  temperature  greater  than  critical. 
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NOMENCLATURE 

a  Exterior  radius,  half-thickness  of  slab,  cm 
c  Specific  heat  capacity,  cal/g-°C 
d  Differential  operator 

E  Activation  energy,  cal/mole  talso  in  kcal/mole  on  nomographs) 

F  Function  Xt^/pca^ 

F(  )  One  function  of  the  quantity  in  parenthesis 

/(  )  One  function  of  the  quantity  in  parenthesis 

h  Heat-transfer  coefficient,  cal/cm^-sec-°C 
k  Parameter  in  Eq.  3 

log  Logarithm  to  base  10 

V  Modulus  cm,  unit  of  linearly  plotted  variable 
QZ  Heat  release  parameter,  cal,'  g-sec 
R  Universal  gas  constant,  1.987  cal,  mole-  K 
r  Radius,  distance  from  central  plane  of  slab,  cm 
T  Temperature,  °K 

Critical  temperature,  °K 
Tj  Temperature  of  surface  of  object,  °K 
To  Temperature  of  surroundings,  °K 

Tl  Constant  temperature  of  surface,  held  at  value  greater  than  critical  temperature,  °K 
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7*  l/(l/7„  -  l/7i).  °K 

t  Time,  sec 

tg  Time-to-explosion,  sec 

Variables  in  formulas  for  types  of  nomographs 
Z  See  QZ 

a,/3  Axes  on  nomographs 

B  Shape  factor,  Eq.  4.  Numerical  values  herein  are  critical  values  of  5,  see  Eq.  9 
0  Dinensionless  variable  £(7  —  T^/RT^ 

Of  Value  of  0  at  surface 
A  Thermal  conductivity,  cal/cm-sec-°C 
p  Density  gm/cm® 

^(  )  One  function  of  the  quantity  in  parenthesis 
^(  )  One  function  of  the  quantity  in  parenthesis 
(U  Dimensionless  distance,  r/a 
S  Partial  differential  operator 

Laplacian  operator  of  vector  analysis 
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INTRODUCTION 

It  has  been  known  for  some  time  that  many  materials,  including  propellants  and  hi ^  explo¬ 
sives,  can  exhibit  the  phenomenon  of  self-heating.  This  is  the  liberation  of  energy  in  die  interior 
of  the  material  caused  by  a  slow  chemical  reaction  that  is  often  a  decomposition  reaction.  this 
heat  is  not  removed  from  the  material,  the  temperature  of  the  material  increases  with  a  concurrent 
increase  in  the  rate  of  reaction.  If  the  material  is  not  disturbed,  a  state  of  dynamic  equilibrium 
will  be  reached  at  which  heat  is  removed  as  fast  as  it  is  generated;  or,  if  heat  cannot  be  removed 
sufficiently  fast,  the  temperature  increases  more  and  more  rapidly,  and  eventually  the  materia 
catches  fire  or  explodes  if  it  is  capable  of  so  doing. 

The  quantitative  description  of  this  type  of  situation  has  been  relatively  recent.  The  differ¬ 
ential  equation  describing  the  situation  is  well  known  and  may  be  written  as 

pc—  =  k\^T  +  pQZ  exp(-£ A* T)  (1) 

dt 

in  which  it  is  assumed  that  the  thermal  conductivity.  A,  is  independent  of  position  and  that  the 

energy  liberated  by  chemical  reaction  is  given  in  terms  of  cal/cm  -sec  by  the  expression 

pQZ  exp(-E/RT).  This  partial  differential  equation  is  strongly  nonlinear  and  has  not  been  solved 

until  recently. 

If  a  steady-state  temperature  distribution  can  be  attained,  it  will  satisfy  the  ordinary  non¬ 
linear  differential  equation 

+  p(^Z  exp{—E/RT)  =0 

Frank-Kamenetski  (lief.  1)  obtained  the  steady-state  solution  for  a  slab  with  both  faces  held 
at  a  constant  temperature;  Chambre  (Ref.  2)  gave  the  steady-state  solutions  for  the  cylinder  and 
the  sphere  with  their  outer  surfaces  held  at  a  constant  temperature;  and  Enig  (Ref.  3)  presented 
the  steady-state  solution  for  the  hollow  cylinder  with  its  surfaces  held  at  constant  temperatures. 
Longwell  and  Conrad  (Ref.  4)  presented  the  solution  for  the  hollow  cylinder  with  its  extenor  sur¬ 
face  held  at  a  constant  temperature  and  its  interior  surface  adiabatic.  The  above  solutions  were 
all  analytical  and  employ  an  approximation  that  represents  the  exponential  of  the  reciprocal  tem¬ 
perature  as  an  exponential  of  a  direct  temperature  function.  Enig,  Shanks,  and  Southworth  (Ref.  5) 
made  numerical  solutions  for  the  steady  state,  which  did  not  use  this  approximation.  Constant- 
temperature  surfaces  were  taken  for  the  slab,  the  cylinder,  and  the  sphere,  and  their  results 
agreed  well  with  those  found  earlier  for  situations  ordinarily  encountered. 

The  foregoing  authors  all  used  constant-temperature  surfaces.  This  implies  relatively  large 
heat-transfer  coefficients.  Thomas  (Ref.  6)  obtained  solutions  for  die  steady  state  for  the  slab, 
the  cylinder,  and  the  sphere  as  functions  of  the  heat-transfer  coefficients  to  the  surroundings  that 
were  taken  to  be  at  constant  temperature. 

Gray  and  Harper  (Ref.  7I  assumed  that  the  thermal  conductivity  was  large  in  comparison  with 
the  heat-transfer  coefficient,  and  they  neglected  temperature  variations  in  the  solid.  Thus  they 
were  able  to  obtain  both  steady-state  and  transient  solutions  for  the  average  temperature  of  the 
material. 

Recently,  Zinn  and  Mader  (Ref.  8)  attacked  the  full  partial  differential  equation  (Eq.  D  by 
numerical  means,  using  high-speed  computing  equipment.  They  worked  w  ith  the  slab,  the  cylinder. 
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and  the  sphere  with  the  surfaces  held  at  sufficiently  high  constant  temperatures  that  a  steady 
state  was  not  possible,  and  they  determined  the  transient  temperature  distributions  and  the  time 
that  would  elapse  before  the  material  exploded. 

Using  the  work  of  the  authors  cited,  it  is  now  possible  to  determine  the  steady-state  temper¬ 
ature  distributions  and,  more  important,  the  maximum  conditions  for  which  a  steady-state  solution 
is  possible  for  several  combinations  of  geometry  and  boundary  conditions.  In  some  circumstances 
it  becomes  necessary  to  exceed  the  maximum  conditions  for  which  a  steady-state  solution  is  pos¬ 
sible,  and,  thus,  if  these  conditions  were  exceeded  for  a  sufficient  length  of  time,  the  material 
would  explode.  It  is  now  possible,  for  constant-temperature  surfaces  on  the  slab,  the  cylinder, 
and  the  sphere  (corresponding  to  very  large  heat-transfer  coefficients),  to  determine  how  long  the 
material  will  take  to  explode.  This  information  is  important  to  engineers  who  process  explosives 
or  design  weapons. 


DESCRIPTION  OF  THE  PROBLEM 


The  differential  equation  for  the  steady  state,  Eq.  2,  can  be  put  into  dimensionless  form  as 


d^d  lic\/dd\ 

"  (d(  J  ^ 

(3) 

in  which 

pQZEa^ 

(4) 

It 

1 

o 

(5) 

(I)  =  r/ a 

(6) 

and  A  =  0,  1,  and  2,  respectively,  for  the  slab,  the  cylinder,  and  the  sphere, 
of  the  type 

Boundary  conditions 

h(T,  -  To)  +  \{dT/dr),  =  0 

(7) 

become 

(Aa/A)(?,  +  (dd/dcj),  =  0 

(8) 

If  the  dimensionless  ratio  Aa/A  is  zero,  the  surface  is  adiabatic;  if  the  ratio  is  infinite,  the 
surface  assumes  the  temperature  of  the  surroundings,  T q.  For  intermediate  ratios,  the  surface 
temperature  is  somewhat  higher  than  that  of  the  surroundings. 

The  critical  temperature  of  an  object  is  defined  as  the  maximum  temperature  of  the  surround¬ 
ings  imposed  on  the  object  in  some  specified  fashion  for  which  a  steady  state  is  poss  lde.  It  fol¬ 
lows,  therefore,  that,  if  the  surroundings  are  held  at  a  temperature  greater  than  the  critical  temper¬ 
ature  for  the  specified  situation,  the  object  will  self-heat  until  it  ignites  or  explodes. 

The  maximum  or  critical  values  of  the  parameter  S  have  been  given  by  several  authors  (Ref. 
1—6)  for  various  geometries  and  boundary  conditions.  Values  for  the  slab,  the  solid  cylinder,  and 
the  sphere  have  been  given  by  Frank-Kamenetski  and  Chambre.  These  correspond  to  boundary 
conditions  with  infinite  heat-transfer  coefficients  and  represent  values  of  S  of  0.88,  2.00,  and 
3.32,  respectively. 

Table  1  shows  values  of  8  for  hollow  cylinders  as  functions  of  the  ratio  of  inner-  to  outer- 
radius.  Two  sets  of  boundary  conditions  are  included.  The  first  is  for  the  same  surface  tempera- 
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TABLE  1.  Shape  Factors  FOR  THE  Hollow  Cylinder 


Shape  factor,  8 

Radius 

ratio 

Same  temperature 
inside  and  outside 

Adiabatic  inside, 
constant  temperature  outside 

0.025 

3.008 

2.003 

0.05 

3.320 

2.020 

0.10 

3.892 

2.077 

0.15 

4.492 

2.172 

0.20 

5.171 

2.306 

0.25 

5.966 

2.483 

0.30 

6.921 

2.713 

0.35 

8.091 

3.004 

0.40 

9.556 

3.384 

0.45 

11.43 

3.878 

0.50 

13.88 

4.530 

0.55 

17.19 

5.413 

0.60 

21.81 

6.646 

0.65 

28.55 

8.434 

0.70 

38.91 

11.17 

0.75 

56.10 

15.69 

0.80 

87.73 

23.92 

0.85 

156.1 

41.57 

0.90 

351.3 

91.50 

0.95 

1405 

358.4 

1.00 

1 

J - - - 

eo 

ture  inside  and  out;  the  second  is  adiabatic  inside  and  constant  surface  temperature  outside. 

Thus,  the  first  set  of  boundary  conditions  corresponds  to  infinite  heat-transfer  coefficients  an 
the  same  temperature  of  surroundings  inside  and  out.  The  values  of  B  are  those  given  by  Enig  in 
Ref.  3.  The  second  set  of  boundary  conditions  corresponds  to  an  infinite  heat-transfer  coefficient 
on  the  outside  and  a  zero  heat-transfer  coefficient  on  the  inside.  The  values  of  8  for  radius  ratios 
less  than  0.35  were  calculated  using  equations  given  by  Conrad  and  Longwell  in  Ref.  4,  and 
those  for  the  remaining  ratios  were  derived  from  the  results  of  Enig  as  described  previously. 

Table  2  gives  values  of  S  for  finite  heat-transfer  coefficients  for  the  slab,  cylinder,  and 
sphere.  These  values  were  read  from  graphs  presented  by  Thomas  in  Ref.  6  and  are  given  as 
functions  of  the  dimensionless  ratio  ha/K 

Knowing  S,  the  critical  temperature  may  be  calculated  from  Eq.  4,  which  can  be  arranged  to 

give 

_ E  _  (9) 

"  2.303«  log  [pa'^QZE/kRTlS) 

The  temperatures  are  absolute  temperatures  expressed  in  degrees  Kelvin,  and  the  activation 
energy,  E,  is  in  cal/mole;  not,  as  it  is  usually  quoted,  in  kcal/mole.  The  units  and  definitions  of 
the  various  quantities  are  listed  in  the  nomenclature. 

Zinn  and  Mader  (Ref.  8)  presented  a  graph  showing  the  dimensionless  parameter  f  =  Xl^/pca 
as  a  function  of  £(1/ -  \.  Tx)  for  the  slab,  the  cylinder,  and  the  sphere.  The  boundary  condi¬ 
tions  are  that  the  surface  is  held  at  the  temperature  Ij,  which  is  above  the  critical  temperature. 
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TABLE  2.  Shape  Factors  With  Finite 
Heat-Transfeh  Coefficients 


Dimensionless 

ratio, 

ha/\ 

Shape  factor, 

s 

Slab 

Cylinder 

Sphere 

0 

0 

0 

0 

1 

0.26 

0.55 

0.87 

2 

0.42 

0.91 

1.47 

3 

0.51 

1.13 

1.86 

4 

0.57 

1.28 

2.13 

5 

0.61 

1.39 

2.31 

6 

0.63 

1.47 

2.43 

8 

0.68 

1.58 

2.59 

10 

0.71 

1.67 

2.74 

12 

0.74 

1.73 

2.83 

14 

0.76 

1.77 

2.91 

OC 

0.88 

2.00 

3.32 

and  thus  correspond  to  an  infinite  heat-transfer  coefficient  to  surroundings  at  Ti.  In  order  to  cal¬ 
culate  a  time-to-explosion,  t*,  the  critical  temperature  is  first  calculated  from  Eq.  9,  then  the 
quantity  £(l/7„  —  \./T-\)  is  calculated,  and  F  read  from  the  graph.  Finally,  is  calculated  from 
the  value  of  F. 


USE  OF  NOMOGRAPHS 

The  calculations  of  the  critical  temperature  and  of  the  tinie-to-explosion  can  be  done  easily 
by  the  use  of  nomographs  that  have  been  constructed.  These  nomographs  have  the  usual  advan¬ 
tages,  such  as  rapidity  of  use,  reduction  of  the  likelihood  of  making  errors,  and  the  presentation 
of  results  in  desired  units  and  can  be  used  by  persons  who  are  not  thoroughly  conversant  with  the 
calculations  represented.  The  construction  of  these  nomographs  is  described  in  the  Appendix; 
their  use  is  described  below. 

CRITICAL  TEMPERATURE 

The  critical  temperature  of  an  object  is  the  maximum  temperature  of  the  surroundings  imposed 
on  the  object  in  some  specified  fashion  for  which  a  steady  state  is  possible.  In  order  to  calculate 
the  critical  temperature  using  Nomograph  1,  seen  in  Fig.  1,  four  quantities  are  needed: 

1.  a,  the  half-thickness  of  a  slab,  or  the  exterior  radius  of  a  cylinder  or  sphere, 
expressed  in  cm. 

2.  The  quantity  pQZ,E/\R,  expressed  in°KVcm^.  Note  that  E  must  be  expressed 
in  cal/mole  in  calculating  this  quantity.  The  quantity  QZ  is  in  cal/'g-sec. 

3.  E,  the  activation  energy,  expressed  in  kcal/mole. 

4.  The  critical  value  of  the  shape  factor  5  for  the  specific  geometry  and  boundary 
conditions.  Values  of  S  can  be  obtained  as  follows: 

a.  Slab,  solid  cylinder,  or  sphere,  with  infinite  heat-transfer  coefficient  (constant- 
temperature  surface):  use  points  labeled  “slab,”  “cylinder,”  or  “sphere”  on  the  S 
axis  of  Nomograph  1. 

b.  Hollow  cylinder,  infinite  heat-transfer  coefficients  inside  and  outside,  same 
temperature  inside  and  outside:  use  Table  1,  column  headed  “Same  temperature 
inside  and  outside”  for  proper  ratio  of  inside-to-outside  radii. 


NAVWEPS  REPORT  7646 


NAVWEPS  REPORT  7646 


c.  Hollow  cylinder,  infinite  heat-transfer  coefficient  outside,  zero  heat-transfer 
coefficient  inside:  use  Table  1,  column  headed  “Adiabatic  inside,  constant 
temperature  outside”  for  proper  ratio  of  inside-to-outside  radii. 

d.  Slab,  solid  cylinder,  or  sphere,  with  finite  heat-transfer  coefficient:  calculate 
ha/X,  take  value  of  S  from  Table  2. 

The  use  of  Nomograph  1  is  shown  by  the  key.  A  straight  line  drawn  between  a  and  pQZE/XR 
intersects  the  a  axis  at  a  point  a,  a  straight  line  between  the  point  ol  and  the  value  of  8  intersects 
the  axis  at  a  point  and  a  straight  line  between  the  point  and  the  value  of  E  intersects  the 
axis  at  a  point  that  is  read  as  the  critical  temperature,  T^.  On  the  nomographs,  the  temperatures 
are  expressed  in  degrees  centigrade  in  order  to  facilitate  calculations. 

This  will  be  demonstrated  by  four  examples,  one  for  each  of  the  four  general  types  of  situa¬ 
tions  presented  above  for  shape  factors  S. 

Example  1,  A  solid  cylinder  20  cm  in  diameter,  composed  of  material  having  £  =  55  kcal/ 
mole  and  pQZE/XR  =  1.0  x  10*°  °K*/cm^,  has  its  exterior  surface  held  at  constant  temperature. 
Presume  the  length-to-diameter  ratio  is  large,  so  the  end  effects  are  negligible,  and  the  point 
labeled  “cylinder”  on  the  S  axis,  corresponding  to  S  «=  2.00,  is  appropriate.  Using  a  =  10  cm  and 
the  other  data  as  given.  Nomograph  1  gives  a  critical  temperature,  Tr.,  of  182°C.  The  construction 
is  shown  dotted  on  the  nomograph.  A  numerical  calculation  using  Eq.  9  gives  455. 6°K  (182.4‘T;). 

Example  2.  A  hollow  cylinder  6.0  inches  in  outside  diameter  and  3.9  inches  in  inside  diam¬ 
eter  is  held  at  the  same  constant  temperature  on  both  surfaces.  E  =  44  kcal/mole  and  pQZE/XR 
=  5  X  10*^  °K*/cin*.  Presume  the  length-to-dianieter  ratio  is  large,  so  Table  1  can  be  used.  The 
radius  ratio  is  3. 9/6.0  =  0.65  and,  from  Table  1,  S  =  28.55.  a  =  6.0  x  2.54/2  •=  7.62  cm.  The  crit¬ 
ical  temperature  is  read  from  Nomograph  1  as  149°C.  As  a  check,  the  numerical  calculation  using 
Eq.  9  gives  422.5°K  (149.3°C). 

Example  3.  The  hollow  cylinder  of  Example  2  is  held  at  a  constant  temperature  on  the  out¬ 
side,  but  there  is  no  heat  transfer  on  the  inside.  Table  1  is  used  with  the  radius  ratio  of  0.65, 
giving  8  5=  8.434.  a  =  7.62  cm,  £  =  44  kcal/mole,  and  pQZE/XR  =  5  x  10*^°K*/cm^,  The  nomo¬ 
graph  gives  139°C  as  the  critical  temperature.  As  a  check,  the  numerical  calculation  gives 
412.6°K  (139.4°C). 

Example  4,  A  sphere  20  cm  in  diameter  is  suspended  in  an  air  stream,  and  the  heat-transfer 
coefficient  is  estimated  at  2.0  Htu/ft*-hr-°F.  The  thermal  conductivity  of  the  material  is  0.2 
Btu/hr-ft-°F.  £  =  55  kcal/mole,  and  pQZE/XR  =  1.0  x  10*°  °K*/cm*.  a  is  10  cm.  The  parameter 
ha/X  must  be  calculated,  watching  units  so  that  the  result  is  dimensionless.  In  this  case,  h  must 
be  converted  to  cal/cm*-sec-°C  and  k  to  cal/cm-sec-°C,  or  the  radius  converted  to  feet.  Doing 
the  latter,  a  —  10/(12  x  2.54)  =  0.328  foot  and  ha/X  =  2.0(0. 328)/0. 20  =  3.28.  Using  Table  2,  8  is 
interpolated  as  1.94.  The  nomograph  gives  182°C.  The  numerical  calculation  gives  455.3°K 
(182.1°C). 

TIME-TO-EXPLOSION 

The  time-to-explosion,  can  currently  only  be  calculated  for  the  slab,  the  solid  cylinder, 
and  the  sphere  and  with  the  outside  surfaces  held  at  a  constant  temperature  corresponding  to  an 
infinite  heat-transfer  coefficient.  This  constant  temperature,  designated  as  Zj,  is  larger  than  the 
critical  temperature  for  the  infinite  heat-transfer  coefficient.  The  time-to-explosion,  is  calcu¬ 
lated  using  Nomographs  2  and  3  (Fig.  2  and  3),  and  the  following  data  nre  required  in  addition  to 
those  required  to  calculate  the  critical  temperature  using  Nomograph  1. 

1.  1 1,  the  constant-temperature  surface,  °C. 

2.  'ITie  quantity  pc/X,  scc/ cni^. 
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The  critical  temperature,  T„,  is  calculated  for  the  constant-temperature  surface  as  described 
above  and  shown  in  Example  1.  Nomograph  2  is  used  as  shown  by  the  key.  It  is  entered  with 
and  Jj,  which  are  connected  by  a  straight  line  intersecting  the  o.  axis  at  a  point  a.  The  point  a 
is  connected  with  £  by  a  straight  line  that  intersects  the  E/T*  axis  at  a  point  where  the  value  is 
read  as  E/T*.  This  value  is  taken  to  Nomograph  3,  and  the  corresponding  point  on  the  appropriate 
E/T*  axis  (sphere,  cylinder,  or  slab)  is  connected  to  the  aaxis  by  a  straight  line  through  the 
point  marked  ®.  The  intersection  on  the  ot  axis  is  connected  with  the  value  of  a  on  the  a  axis  to 
to  give  an  intersection  on  the  /3  axis,  and  the  latter  is  connected  with  the  value  oi  pc/\  by  a 
straight  line  that  jntersects  the  t^axis  at  the  time-to-explosion. 

This  calculation  is  illustrated  by  an  example. 

Example  5.  A  long,  solid  cylinder  10  cm  in  radius,  composed  of  materials  having  £  =  55 
kcal/mole,  pQZE/XR  =  1.0  x  lO^^^KVcm^,  and  pc/X  =500  sec/cm^  has  its  exterior  surface 
held  at  240°C.  The  critical  temperature  for  the  constant-temperature  surface  was  found  in  Exam¬ 
ple  1  to  be  182°C.  Using  T^  =  182°C  and  Tj  =  240°C,  the  value  of  E/T*  is  read  from  Nomo¬ 
graph  2  as  13.7  cal/mole-°K.  This  value  is  used  on  the  cylinder  E/T*  axis  of  Nomograph  3  and, 
using  the  given  values  of  a  and  pc/X.  the  time-to-explosion  is  read  as  2,500  seconds.  As  a  check 
by  numerical  calculation,  T^  was  found  to  be  182. 4°C  in  Example  1.  Using  this  value  and  the 
given  values  of  Tj  and  £,  E/T*  is  calculated  as  13.55  cal/mole-°K.  Again  using  this  calculated 
value  and  the  given  values  of  a  and  pc/X,  tg  is  found  to  be  2,450  seconds. 


COMMENTS  AND  CAUTIONS 

The  nomographs  accompanying  this  report  allow  fairly  accurate  solutions  of  certain  mathe¬ 
matical  problems,  and,  if  physical  problems  are  well  described  by  the  mathematics  involved,  the 
solutions  of  the  physical  problems  can  be  equally  well  determined.  As  a  matter  of  fact,  the  phys¬ 
ical  problems  are  not  always  as  simple  as  the  mathematical  problem  described  by  Eq.  1  and  the 
associated  initial  and  boundary  conditions. 

The  major  trouble  is  that  the  reaction  rate  may  not  be  well  described  by  an  expression  of  the 
type  QZ  exp(—E/KT),  in  which  all  parameters  except  the  temperature  are  constants.  Decomposi¬ 
tion  reactions  are  sometimes  autocatalytic  or  catalyzed  by  small  amounts  of  impurities,  so  that 
the  decomposition  rate  is  a  function  of  duration  or  extent  of  decomposition,  and,  if  this  is  true  to 
a  significant  extent,  the  analysis  described  herein  is  not  suitable.  If  an  autocatalytic  reaction, 
suck  as  the  decomposition  of  nitrocellulose,  is  controlled  by  the  use  of  stabilizers,  the  relations 
used  can  only  be  valid  as  long  as  the  stabilizer  remains  and  is  effective. 

Appropriate  values  of  the  physical  constants  must  be  known.  Determinations  of  thermal  con¬ 
ductivity,  specific  heat,  and  density  can  be  relatively  precise,  but  the  measurement  of  the  reaction- 
rate  parameters  QZ  and  £  is  difficult,  and  the  accuracy  of  such  parameters  is  generally  much 
poorer.  Generally,  the  accuracies  of  these  two  parameters  are  the  limiting  factors  on  the  accuracy 
of  determination  of  critical  temperatures  or  times-to-explosion,  even  though  the  decomposition  re¬ 
actions  do  appear  to  be  well  behaved. 

No  phase  changes  have  been  considered  in  this  treatment.  If  phase  changes  occur,  the  criti¬ 
cal  temperature  should  not  be  affected  appreciably,  but  the  time-to-explosion  will  probably  be 
somewhat  greater  than  calculated. 

As  a  result,  the  values  of  critical  temperature  and  time-to-explosion  determined  should  be 
considered  as  only  approximate.  They  should  be  used  with  margins  of  safety  dictated  by  engi¬ 
neering  judgments  that  are  strongly  influenced  by  the  magnitude  of  the  cost  of  being  wrong. 
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Appendix 

DETAILS  OF  CONSTRUCTION  OF  NOMOGRAPHS 


NOMOGRAPH  1 

This  nomograph  solves  the  equation  for  the  critical  temperature. 

E 


2.303R  logipa^QZE/kRT^S) 


Rearranging  Cq.  10  gives 


303/? 


which  can  be  put  in  the  form  of  Type  C  of  Davis  (Ref.  9)  as 

f(.x)  =  tp(x)  ‘  F  (y)  +  (^(z) 

with 

f(x)  =  27„  log7„ 

2 

F  ' 


^Pix)  =  7„ 


(10) 


(11) 


(12) 


(13) 


2.303/? 


The  value  of  F (y)  is  calculated  using  logarithmic  charts  described  by  Davis  as  Type  B, 
which  are  really  addition  charts  (Davis  Type  A)  in  terms  of  logarithms.  The  modulus  of  a  partic¬ 
ular  scale  is  defined  as  the  distance  along  the  scale  for  one  unit  of  the  linearly  plotted  variable. 
For  example,  in  Nomograph  1,  the  right-hand  scale,  Icbeled  a,  is  actually  a  linear  scale  in  terms 
of  log  and  the  modulus  of  this  scale,  Wa,wa8  taken  as  10  cm  per  unit  change  in  log  a*. 

Referring  to  Fig.  4,  which  shows  an  addition  chart  solving  the  equation, 


fix)  =  Hy)  +  <^(r) 


(14) 


CT 


FICs  4.  Addition  Chari. 
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it  is  necessary  that 


My+M^ 


(15) 


and  that 


AB  My 

Jc 


(16) 


Thus,  the  selection  of  two  moduli  automatically  establishes  the  modulus  and  relative  posi 
tion  of  the  third  scale. 


Nomograph  1  proceeds  as  fdllows: 

First,  <x  is  calculated  as 

>17' 

with  Ma  =  10  cm  and  the  modulus  of  (.pQZE)/{\R)  «  2.5  cm,  giving  Ma.  =  2.0  cm  and  scale  spacing 
from  a  in  the  ratio  of  4  to  1,  which  was  used  as  12  cm  and  3  cm. 

Next,  /3  w.is  calculated  as 

/3  =  a  +  log  (1/5)  =  a  -  log  5  (18) 

with  Mo.  =  2.0,  Ms  =  10.0,  giving  M^  =  1.667  and  scale  spacing  from  /3  in  the  ratio  of  1  to  5,  which 
was  used  as  6  cm  and  30  cm.  Thus, 


^  =  log 


Fiy) 


(19) 


The  recurrent-variable  equation  shown  above  as  Eq.  11  and  12  may  now  be  solved.  The  con¬ 
struction  is  shown  in  Fig.  5.  The  point  F,  which  is  located  on  the  line  connecting  the  given  val¬ 
ues  of  ^(z)  and  F{y),  is  the  location  for  the  desired  value  of  /(*).  It  is  located  in  accordance 
with  the  relations 


CE 

CA 


M, 

Mf  l/*  (*)  +  My 


(20) 


FIG.  5.  Recurrent-Variable  Chart. 
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EF 


[ — - \{h 


■) 


(21) 


Points  C  and  A  are  located  where  ^(z)  and  F  (y)  are  respectively  equal  to  zero,  and,  thus, 
the  line  CEA  is  between  the  zeros  of  the  two  scales. 


For  Nomograph  1,  My  =M^  =  1.667,  and  M^  =  4.575  x  10  where  M^  is  in  terms  of  E /2,303R, 
The  scale  on  Nomograph  1  is  marked  in  terms  of  E,  using  the  modulus  of  1.00  x  10  ®  cm/kcal- 
mole,  which  is  consistent  with  M j.  The  spacing  between  the  E  scale  and  the  /3  scale  was  made 
24  cm,  and,  taking  the  origin  of  Cartesian  coordinates  (x,y)  at  the  60  kcal/mole  point  on  the  E 
scale  and  with  the  other  scales  as  shown  on  Nomograph  1,  the  points  for  the  critical  temperature, 
T„  (expressed  in  degrees  Kelvin),  are  located  at  (x,y)  where 


247-„ 

+  364.4 


(22) 


and 

y  =  60  -  3.8611X:  +  0.13888*  log  T„  (23) 

The  scale  was  graduated  in  terms  of  degrees  Centigrade  and  the  E  scale  in  terms  of  kcal/ 
mole,  although  these  variables  must  be  expressed  in  terms  of  degrees  Kelvin  and  cal/mole,  re¬ 
spectively,  in  the  equations. 


NOMOGRAPH  2 


This  nomograph  calculates  E/T*  from  £,  the  critical  temperature  T„  obtained  from  Nomograph 
1,  and  the  constant-temperature  surface  7;,  according  to  the  equation 


£/7*  =£(l/7„-l/7i) 

(24) 

First,  an  addition  chart  (Davis  Type  A)  is  used  to  calculate 

l/7*=(l/7„-l/7i) 

(25) 

which  is  actually  done  in  the  form 

l/7x  =  l/7„-a 

(26) 

Moduli  satisfying  Eq.  15  were  selected  as  Mt„  ”  2.0  x  10  *  cm-°K,  4.667  x  10*  cm-°K,  and 
Mti  =  1"4  X  10*  cm-°K,  and  the  distances  between  scales  7^  to  7i  and  7i  to  a  as  7.5  cm  and 
17.5  cm,  respectively.  The  scales  T„  and  Ti  were  calculated  as  inverse  absolute  temperatures 
and  graduated  in  degrees  Centigrade.  The  scale,  a,  calculated  as  inverse  absolute  temperature, 
increases  in  the  direction  opposite  to  the  other  two  scales. 

Next,  E/T*  is  calculated  from 

E/T*  =  {a){E)  (27) 

The  scale  for  a  is  linear,  not  logarithmic,  so  a  nonlogarithmic  multiplication  chart  of  Type  D 
as  described  by  Davis  was  employed.  This  type  is  constructed  as  shown  in  Fig.  6  for  the  equation 

f(x)  =  F{y)/^{z)  (28) 

The  scales  for  £(y)  and  <f>{z)  are  laid  out  in  opposite  directions  using  the  moduli  My  and  Af,, 
starting  with  zeros  at  points  A  and  C,  respectively.  The  x  azis  connects  points  A  and  C.  A  point, 
H,  is  selected  on  the  z  axis  with  CH  =  h,  and  a  temporary  auxiliary  scale  with  modulus  M  =  hMy/M, 
is  constructed  on  the  y  axis  starting  at  A.  The  x  scale  on  the  *  axis  is  constructed  by  projecting 
the  auxiliary  scale  by  means  of  straight  lines  through  the  point  H. 

A  disadvantage  of  this  type  of  nomograph  is  that  the  scale  for  /(*)  is  likely  to  be  quite  com¬ 
pressed.  Several  permutations  of  Eq.  27  are  possible,  and  it  was  found  that  the  largest  scales 
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FIG.  6.  Nonlogarithmic  Multiplication  Chart. 


could  be  obtained  by  using  /(*)  =  E/T*,  F(y)  •=  «,  and  =  l/E.  Ratber  than  actually  construct¬ 
ing  the  auxiliary  scale  and  doing  the  construction  by  projection,  an  algebraic  expression  for  the 
location  of  points  on  the  *  scale  was  derived.  Taking  the  origin  at  point  A  of  Fig.  3,  letting  y  now 
be  the  vertical  Cartesian  coordinate  of  the  point  f(x)  on  the  *  axis,  and  with  the  elevation  of  point 
C  above  the  horizontal  axis  being  CJ,  the  expression  is 

Z  _  .  A*)  (29) 

CJ  +  MyF(y)  f(x)  +  M^/My 

Taking  My  =  A/a  '  4.667  x  10*  cin-°K,  and  M,  =  1.8  x  10®  cm-kcal/mole,  and  CJ  =  60  cm,  Eq. 

29  becomes  * 

60(F^  (30) 

^  (E/7*)  +  38.57 

which  was  used  to  compute  the  E/T*  scale  of  Nomograph  2.  The  E  scale  was  calculated  in  terms 
of  l/£  in  mole/cal,  but  graduated  as  E  in  kcal/mole. 

NOMOGRAPH  3 

Zinn  and  Mader  presented  a  graph  of  the  quantity  F  -  XTg/pea^  as  a  function  of  the  quantity 
E/T*,  showing  curves  for  the  slab,  the  cylinder,  and  the  sphere.  Nomograph  3  is  entered  with 
E/T*  and  determines  the  value  of  logF  as  a,  then  calculates  t,  from  the  equation 

1.  -  F(a^)(pc/X)  (31) 

The  first  operation  amounts  essentially  to  reading  the  graph  to  determine  F  for  the  appropriate 
geometry  and  given  value  of  E/T*. 

The  a  scale  was  assigned  a  modulus  (A/a)  of  8.0  cm.  A  swing  point  was  located  midway  be¬ 
tween  the  slab  axis  and  the  cx  axis,  10  cm  from  each.  The  slab  axis  thus  has  a  modulus  of  8.0  cm. 
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The  cylinder  axis  was  located  3  cm  to  the  left  of  the  slab  axis  and  has  a  modulus  of  8.0(10  +  3)/ 
10  =  10.4;  the  sphere  axis,  located  3  cm  to  the  left  of  the  cylinder  axis,  has  a  modulus  of  8.0 
(10  +  6)/10  =  12.8.  In  order  to  make  a  increase  upward,  the  scales  for  the  three  shape  axes  in¬ 
crease  downward  in  terms  of  a;  however,  they  are  graduated  in  terms  of  E/T*,  which  increases 
upward. 

Measuring  y  upward  from  a  horizontal  base  line  intersecting  the  a  axis  at  —3.00,  the  swing 
point  was  located  at  y  =  14.0  cm,  and  the  scales  on  the  three  shape  axes  then  were  graduated  in 
terms  of  E/T*  at  the  locations 

y  =  28.0-8.0  (logf  +  3.0)  slab 

y  =  32.2  —  10.4  (log  F  +  3.0)  cylinder  (32) 

y  =  36.4  -  12.8  (log  F  +  3.0)  sphere 

in  which  the  values  of  logf  were  read  from  the  appropriate  curve  on  a  large  scale  graph  supplied 
by  Mader.  ^ 

The  next  operation  is  to  calculate  /3  from 

/3  =  logfa^  =  a  +  loga^  (33) 

This  is  done  by  addition.  The  a  scale  has  a  modulus  of  10.0  cm  in  terms  of  loga^,  and,  as  Ma.- 
8.0  cm,  is  4.444  cm,  and  the  spacings  of  the  a  and  a  axes  from  the  /3  axis  are  4  and  5  cm,  re¬ 
spectively.  The  a  scale  was  graduated  in  terms  of  a  at  the  positions  given  by  the  equation 

y  =  20  log  a  —  6  (34) 

The  a  scale  and  the  /3  scale  are  not  graduated,  as  these  intermediate  results  need  not  be  read. 
The  final  operation  is  to  calculate  tg  from 

logt*  =  /3  +  log  (pc/X)  (35) 

which  is  also  done  by  addition.  A  modulus  of  the  log  (pc/k)  scale  was  selected  as  5.0  cm,  and, 
as  =  4.444  cm,  Mf  =  2.353  cm  in  terms  of  logt^.  The  /3  and  (pc/k)  scales  are,  respectively,  8 
and  9  cm  from  the  t*  scale.  The  log  (pc/k)  scale  was  graduated  in  terms  of  (pc/A)  at  the  position 

y  =  5  1og(pc/X)  (36) 

The  log  tg  scale  is  appropriately  graduated  in  terms  of  tg  expressed  in  seconds,  and  the  cor¬ 
responding  positions  are  given  by 

y  =  5.647  +  2.353  logtg  (tg,  seconds)  (37) 

However,  it  was  believed  appropriate  to  convert  the  large  values  of  tg  to  larger  units  of  time, 
so  Cq.  37  was  used  for  a  range  of  10~^  to  10*  seconds,  and  the  alternate  equations  following 
were  used  for  a  larger  tg. 


y  =  11.015  +  2.353  log  tg 

(tg,  hours) 

(38) 

y  =  14.263  +  2.353  log  tg 

(tg.days) 

(39) 

y  =  20.292  +  2.353  log  tg 

(tg,  years) 

(40) 

^  Private  communication  to  D.  Lind  from  Charles  L.  Mader,  September  I960. 
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